MOLECULAR RANDOM WALKS AND 
INVARIANCE GROUP OF THE BOGOLYUBOV EQUATION 

O 

S ' YURIY E. KUZOVLEV 

^ . Abstract. Statistics of molecular random walks in a fluid is considered with the 

' help of the Bogolyubov equation for generating functional of distribution functions. 

An invariance group of solutions to this equation as functions of the fluid density is 
discovered. It results in many exact relations between probability distribution of the 
path of a test particle and its irreducible correlations with the fluid. As the conse- 
,-G! ' quence, significant restrictions do arise on possible shapes of the path distribution. 

In particular, the hypothetical Gaussian form of its long-range asymptotic proves to 
be forbidden (even in the Boltzmann-Grad limit). Instead, a diffusive asymptotic is 
allowed which possesses power-law long tail (cut off by ballistic flight length). 



o 



1. Introduction 

Random wandering of particles of the matter is mechanism of diffusion and many 
other transport processes as well as the source of various noises and fluctuations. What 
kind of statistics can rule over it? This important question never was considered in the 
. framework of rigorous statistical mechanics. It may seem that anyway the answer is 



obvious: if in the Lorentz gas statistics of random walks is asymptotically Gaussian pQ 



q 
o. 

Q\ . then all the more it must be Gaussian in usual gas. However, such a reasoning loses its 

O ■ 

convincingness when one notices that it indirectly equates (i) chaotic nature of random 
walk in the sense of Hamiltonian dynamics and (ii) randomness of the walk in the 
sense of the probability theory. The first, as is well known, is characterized by the 
"mixing" of the system's phase trajectories [2], whereas the second by the "statistical 
independence" or simply "independence" [3] of their constituent parts. But already 
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Krylov [2] thoroughly explained that generally the first does not imply the second. 
Therefore the question remains in force. 

It can be comprehended by the example of N hard balls in a box or on torus if 
treating their motion as the motion of single ball in 3iV-dimensional billiard [I] which 
is formed by N(N — l)/2 convex scatterers and resembles the Lorentz gas. How much 
far in time should be apart different fragments of the trajectory of this ball in order 
to behave as statistically independent one on another? It appears that not less than 
by the time necessary for the ball to know about existence of all the scatterers and 
thus about organization of the billiard. Let us express this characteristic time through 
the mean free-flight time in the original three-dimensional system, r ~ (ttTqUVq)^ 1 , 
where tq is diameter of the balls, vq is their typical velocity, v = N/Q their mean 
concentration, and Q is volume of the original system. Since the collision rate in the 
3N -dimensional billiard is equal to summary rate of collisions between arbitrary balls in 
the three-dimensional system, approximately N/2t , then the mentioned characteristic 
time can be estimated as [N(N — l)/2]/[A r /2r] ~ Nt ~ Q/nrQVo . Even for 1 cm 3 of 
the air that is time on order of 1000 years! 

If so, then ergodicity of behavior of iV-particle system can be expected not sooner 
than after observing it during time intervals ^> Nt . From the point of view of real 
physical many-particle systems (let even closed ones), that is quite inaccessible time. 
As to reality, only more or less peculiar parts of the system's phase trajectory are 
observable. Thus, it was true remark [5] that the role of ergodicity in physics is strongly 
exaggerated since there the limit iV — > oo precedes the limit t/r — > oo . 

Similar conclusions follow from the fact that the number of initial conditions what 
determine trajectory of a particle of usual gas is not 6 , as in the case of Lorentz gas, 
but 6iV . Correspondingly, the trajectory can display enough its individuality (and be 
recognized) only after N free flights and N collisions, which requires time ~ Nt . All 
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the more, to disclose its ergodic properties and observe how it divides into statistically 
independent parts we need in time scales ^> Nt . 

To speak more specifically, the non-ergodicity of a system of (infinitely) many par- 
ticles does mean that rates of relaxation and fluctuation processes in it are different 
at different phase trajectories, as well as at diverse parts of one and the same trajec- 
tory [5J. For instance, individual diffusivity of a gas particle undergoes low- frequency 
fluctuations [7j which manifest themselves in a substantial difference of asymptotic 
probability distribution of the particle's displacement from the Gaussian distribution. 

The theoretical tools for rigorous investigations of such a statistics were developed by 
N. Bogolyubov [8] and his followers. That are the Bogolyubov-Born-Green-Kirkwood- 
Yvon (BBGKY) hierarchy of equations and the equivalent Bogolyubov's equation for 
generating functional of many-particle distribution functions. But nobody have learned 
how to use them avoiding truncations of the hierarchy on the hypothesis about some 
"obvious" "statistical independence" adopted from the probability theory. Hence, there 
is necessity of methods which can treat the basic equations respectfully. 

The present work is just caused by such the necessity. Here, we will find and describe 
an invariance group of solutions to the Bogolyubov equation, as applied to the problem 
about statistics of random walk of test particle in thermodynamically equilibrium fluid. 
Then we will discuss consequences from exact relations of the group. In particular, 
demonstrate that they imply non-Gaussian character of long-range asymptotic of the 
random walk. Importantly, in case of gas this asymptotic, along with above mentioned 
characteristic time ~ Nt , is insensible to gas density and stays strongly non-Gaussian 
even in the Boltzmann-Grad limit. 



2. Equations of molecular random walk 

Let a box with volume Q contains N ^> 1 identical atoms plus one more test 
particle. Atoms have mass m, coordinates and momenta p., (j = 1,2... N) and 
interact with each other via potential $ a ( r j — r fc) • The test particle has mass M , 
coordinate R , momentum P and interacts with atoms via potential $&(rj — R) . The 
potentials are spherically symmetric and short-range with impenetrable point core. 
Because of interactions the test particle is in chaotic motion, therefore let us name it 
"molecular Brownian particle" (BP). 

We are interested in probability distribution of current position of BP, R(t) , un- 
der condition that at initial time moment t — it was placed at certainly known 
position: R(0) = Ro , while personal positions of atoms all the times are unknown. 
The simplest statistical ensemble what satisfies this requirement is determined by the 
Liouville equation, dD^/dt = [H^ , -Djv] > f° r full normalized distribution function 
of the system, , and the initial condition to it, 

5(R-Ro) e~ HN / T 
DN(t = 0) = fdRfdPf 1 ...f N 5(R-R )e-^ (2 ' 1} 

where is full Hamiltonian of the system (including interactions with the box walls) 

and f,... = j j ... drkdpk ■ Evidently, such the ensemble differs from the Gibbs 

canonic ensemble by the initial BP's localization only. The latter does not prevent us 

to introduce the marginal distribution functions (DF) F n (t) = Q n f n+1 ... f N D N (t) 

and then go to the thermodynamical limit (iV^oo, Q —>■ oo , v — N/Q = const ) 

in exact analogy with [8]. If writing out all the arguments of DF, we have F n (t) = 

F„(t,R,rM,P,p( n )|R ;z/) , where = {n,...r n }, p (n) = {pi,...p n }. As in the 

Bogolyubov's book, all DF are not normalized in respect to the atoms' coordinates. 

Instead, they must obey the conditions of decoupling of inter-particle correlations under 

separation of particles (in essence, that are conditions of existence of thermodynamical 



limit [8], [9]). For our task, with taking into account complete symmetry of DF in 
respect to atoms, these conditions can be written as follows: F n — ► F n _iG m (p n ) at 
r n — > oo , where G m (p) = (27r Tm)~ 3//2 exp (— p 2 /2Tm) is the Maxwell momentum 
distribution of a particle with mass m . The only, non-principal, difference from [8] is 
that numeration of DF begins from zero, so that F Q (t, R, P| R ; v ) describes the state 
of BP, and in respect to BP's coordinate all the DF are normalized. In particular, 
jF dH = Gm(P) ■ The basic Liouville equation induces the BBGKY equations 

dF d f n d f 

-^=[H n ,F n } + u — J S k '(R - r n+1 ) F n+1 + v ^ — \ - r n+1 ) F n+1 

(2.2) 

(n = 0, 1, . . . ) with initial conditions 

F n (t = 0) = 5(R - Ro) Fi e<?) (r (n) |R; v) G M (P) Y[ n J=1 G m ( Pj ) , (2.3) 

where H n is Hamiltonian of subsystem "n atoms plus BP " , $^ 6 (r) = V$ a ,6( r ) , and 
Fn(r^ |R; z/) are usual thermodynamically equilibrium DF for n atoms in presence 
of BP occupying point R . 

In principle, all that will do for finding Fo(t, R, P| Ro ; v ) and thus probability 
distribution of BP's path, AR(£) = R(£) — R , without any additional assumptions. 
At that, since our problem is expressed by the recurrent relations, it can be lighten, as 
usually, by use of generating functions. Such the approach to the BBGKY hierarchy was 
formulated and sampled already by Bogolyubov [5|. Here it will help us to visualize 
some properties of the hierarchy what are hardly seen directly from it. Following 
Bogolyubov, let us combine all DF into generating functional (GF) 

^{t,R,P,^|Ro;4 = F +X)^r /- Fn LNfo'Pi) ( 2 - 4 ) 

n=l U - Jl Jn j=l 

and equations (\2.2h into corresponding "generating equation" for it: 



where C is operator composed of usual and variational derivatives, 



\:){X\ ) — — ■ 



+ 



(2.6) 



I / [l+#Tl)] 



$ 6 (R-ri), 
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+ 



+ ^// [ i +^(^i)][i+^2): 



$a(ri - r 2 



^(^2) 

with Xj = {rj,pj} . This is direct analogue of equation (7.9) from [8]. To make this 
evident, notice that ip(x) = u(x)/v, where u(x) is functional argument used in [8], 
and factor [1 + if)(xi)] [1 + i ) ( x 2)) can be replaced by [tp(xi) + i J (x 2 ) + ip(xi)ip(x2) ] 
due to the identity J ± L [ $ a (ri — r 2 ) , ... ] = . 

Let us formulate the initial condition to equation (12.51) . For that, it is convenient to 
introduce the generating functional of equilibrium DF in the configuration space, 

R;u} = 1 + jr ^ / F^\r^\ R; v) f[ </>{ Tj ) dr 3 , (2.7) 

n=l ■ J J j=i 

where 0(r) is corresponding functional variable. Besides, introduce linear mapping 
(fi{ip} of functions ip(r, p) to functions 0(r) as follows: 



0{^}(r) 



V»(r,p) G m (p) dp 



{21 



Then the initial conditions (12.31) do intend 

F{0, R,P,^|R ;^} = 5(R-R O )G M (P)^{0W|R;4 



(2.9) 



It is easy to guess or verify that expression G M (P) F {eq) {(j>{ilj}\ R; v] represents the 
stationary solution of equation (I2.5p . so that 



P/M) • d/dR + £ G M (P) F {eq) {<f){^}\ R; u} = 



This equality yields 



d_ 
dr 



~ R) 

T 



dr' , (2.10) 



7 

where, similarly to (|2.7p . <p(r) appears as independent functional variable in the con- 
figuration space. The latter equation determines all the equilibrium DF and is analogue 
of equation (2.14) from [8]. 

Unfortunately, to the best of my knowledge, non-stationary solutions to equations 
like (12.51) above or (7.9) in [8] never were investigated by Bogolyubov or others. How- 
ever, the experience of work with the BBGKY equations prompts (see e.g. [TD]) the 
desirability of a change of variables, i.e. transition from DF to another functions which 
help to concentrate on inter-particle correlations and statistical dependencies. To make 
a suitable choice of such functions, let us discuss hypothetical equalities 

F n (t, R, r<"\ P, pW| R ; v) = F Q (t, R, P| R ; v) F!f q \r^\ R; u) Y\ n j=l G m (Pj) 

or, equivalents, T{t, R, P, | R ; u} = F (t, R, P| R ; v) {(f){i>} \R;u}. They 
state, evidently, that the conditional DF (conditional probability distributions) of 
atoms, i.e. F n /F , do not depend on the past BP's displacement R — Ro . This 
assumption seems natural as concerns thermodynamically equilibrium random walk 
when all possible positions of BP are equivalent. Nevertheless, this is certainly wrong 
assumption, since it is incompatible with the equations 02.21) . Indeed, the substitution 
of Fi in the first equation (for F ) by F (t, R, P| Ro; v) (ri| R; u)G m (pi) turns 
the integral of BP's interaction with atoms into zero: L $ fe '(R — r x ) F\ = (at least 
because F^ (ri|R;i/) is even function of the difference ri — R while $ 6 '(R — ri) is 
its odd function). As the result, the first of equations (12. 2p reduces to the equation of 
free BP's flight, dFo/dt = [Ho , Fq] , as if BP does not interact with atoms at all. 

The aforesaid implies that, first, the current state of the medium (system of atoms) is 
statistically dependent on the summary displacement of BP AR(£) = R(£)— R during 
all time of its observation. In other words, some essential and quantitatively significant 
correlations between the medium and the history of BP's wandering are in existence. 



We will name them "historical correlations" , in order to ideally separate them from the 
equilibrium correlations, as described by F n eg ^ (r^ | R; v) , between current coordinates 
of BP and atoms. Second, to perform an adequate formal separation of the two sorts 
of BP-atoms correlations, we can describe the historical correlations with functions 
V n = V n {t, R, r , P, P^l Ro! v) defined by the generating relations as follow: 

^{t,R,P, ^|R ;4 = V{t,R,P, V;|Ro;4 ^ (e9) {0M|R;*4 , (2.11) 
V{t,R,P,^|Ro;4 = V +J2^7 / ••• / VnJl^Pj) 

n=l U - Jl Jn j=l 

In particular, Vo(t, R, P| Ro; v) = F (t, R, P| R ; v) and 

F 1 (t,R,r 1 ,P,p 1 |R ;i/) = 

= F (t,R,P|R ;z/)F 1 (ef ' ) (r 1 |R;z/)G' m (p 1 ) + Vi(t, R, r x , P, p x | R ; v) 1 ' ' 

It is clear from this definition that from the viewpoint of the probability theory the 
functions V n ( n > ) represent a kind of cumulant functions (cumulants, or semi- 
invariants). Therefore we will name them "cumulant functions" . Notice that in physical 
literature (see e.g. [10]) similar objects frequently are called "correlation functions". 

In terms of the cumulant functions (CF) initial conditions (12.31) and (12.91) become 
strongly simplified: 

V o (0,R,P|R ;z/) = 5(R-Ro)G M (P) , 
K(0,R,rM,P,pW|R ;z/) =0 (n > 0) , 
V{0, R, P, V| R ; u} = 6{R - R ) G M (P) (2.13) 

The conditions of the decoupling of correlations at infinity also become simpler: 

7„(t,R,rW P,p(">|Ro;i/) -> at r k -> oo (2.14) 

( 1 < A; < n ), that is CF tend to zero if at least one n atoms moves away from BP 
to infinity. In opposite, the generating equation (I2.5P in terms of CF becomes more 



complicated. Inserting (I2.1ip into (12.51) one can obtain 
where the new operator appears, 

Z'fa'Jj) = {/[ 1 +^ r )]^( R - r )^{ r ^l R ^} rfr }(]^ + ^) + 
+ J J [l + V(^i)][l + ^a)] $a(ri-r 2 ), z/C{r 2 ,0|R; u }G m (p 2 )— -j 

Here 0(r) and are mentioned as linear functionals of ^(r, p) in the sense of the 
mapping 0{?/>}(r) defined by (12.81) . and besides one more functional is introduced, 

C{r, 1 R; u) = IJ\ 1 (2.16) 

Correspondingly, the BBGKY equations become more complicated. Therefore, here 
we write out them only for extreme but principally interest case of "BP in ideal gas" 
(when atoms do not interact with each other, i.e. $ a ( r ) — ). In this case, 

dV P dV d /• 

77-^7 + ^^77 / ^(R-ri)Vi , 



<9t M dR OP J x 

8Vr 



dV d f 

= [H n ,V n ] + u — J ^' b (R-r n+1 )V n+1 + (2.17) 

n f P 8 \ 

+ T^P(j,n)G m ( Pn )£'(r n -R) - + - K^i (n > 0) , 

where tf n = P 2 /2M+£™ =1 [ pf/2m+$ 6 (R - r,) ] , E(t) = exp [ - $(r)/T] , E'(r) = 
V£?(r) = — $'(r)E(r)/T , and V(j,n) denotes operation of transposition of argu- 
ments Xj = {rj,pj} and = {r n ,p n }. Thus, the natively bidiagonal BBGKY 
hierarchy if being represented via CF becomes al least tridiagonal. 
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3. Invariance group of the equilibrium generating functional 

Further, let us consider the functional (12.161) . which is interesting already because it 
enters the generating evolution equation (12.1 5p . and expand it into a series: 

C{r,0|R;4 = F 1 (e9) (r|R;z/) + ^^ /.../ C n+1 (r, n... r n | R; v) \[ <f>{ rj ) dr,- (3.1) 

n = l n - J J j=i 

According to (12.71) . this expansion just defines the functions C n (ri... r n | R; v) (n > 1 ). 
In particular, 

C 2 (r,n|R ;i /) = F^^r^u) - Ft q \r\K;u) F^ir^u) , 
C 3 (r, ri ,r 2 ) = Ft\r,r u r 2 ) + 2^ ) (r)^(r 1 )F 1 M (r J )- 

_ F M (r ) F N) (ri)r2 ) _ ^(rOF^Cr.ra) - if g) (r 2 ) F 2 (eg) (r, r a ) , 

(for brevity in the second expression arguments R and z/ are omitted). Clearly, from 
the viewpoint of the atoms' coordinates, the functions C n relate to the equilibrium 
DF Fn^ like cumulants of a random field relate to its statistical moments. 

The conditions of vanishing of inter-particle correlations at infinity mean that all 
these cumulants tend to zero if at least one atom is moved away from the others: 

C n {vx ... T n \ R; v) ^ at r., oo (3.2) 



(3.3) 



In case of mowing away BP from atoms we have: 

F[ eq \r\H]v) -> 1 at R -> oo , 
C n (ri ... r n | R; v) -> C n (r 1 ... r n ; v) at R -» oo 

where C n {j\ ... r ra ; v) are cumulant functions of equilibrium media taken in absence 

of BP. Moreover, under sufficiently short-range interaction potentials, all the limit 

transitions in fl3.2H - fl3.3j) are absolutely integrable, which will be assumed below. 

Next, consider equation (12. lOj) rewriting it in the form 

d_ + *6 ; (r-R 



dr T 



C{r, |R; *4 = - J[ 1 + 0(r') ] *'Jj> - r) " } dr' + 

+ C{r, <\> |R; !/} V - j [ 1 + 0(r') ] ^(r' - r) C{r', <\> |R; !/} dr' , (3.4) 
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and with its help let us ascertain (probably, for the first time) some important proper- 
ties of the functional C{r, <fi | R; u} . Notice, first, that integrability of the asymptotic 
(13. 2p makes it possible to extend this functional to bounded functions 0(r) which do 
not tend to zero at infinity, in particular, to constants. This fact allows to introduce 
such the objects as follow: 



00 ii n (j n f f 

lim C{r, 0| R; v} = 1 + — — / ... / C n+1 (r, ri ... r n ; v) , (3.5) 
r)^a ~ n! Ji J n 



R — > 00 



C a {r, 4>\R] v} 



C{r, a + I R; u} 



(3.6) 



where a = const , J ... = J ... dr n , and integrals in (I3.5P are practically independent 
on r . Second, after the change 0(r) — > <r+0(r) the equation (13.41) can be transformed, 
through elementary algebraic manipulations, into equation for C a {r, <p |R; u} : 

5C ff {r',0|R; u} 



d ggr-R) 
9r T 



+ C CT {r,0|R; z/} 



C{r,0|R; u} = - 



T 



1 + 



l + a 



al j 5[0(r)/(l + ( r)] 



1 



l + a 



^(r'-r) C CT {r',0|R;z/} dr' (3.7) 



It is easy to see that its only difference from (13 .4p is the scale transformation of the 
functional argument, 0(r) — > 0(r)/ (l + a) , and besides transformation of the density 
argument v into 



v(a, v) = v C(cr, v) (1 + a) 



(3i 



Third, formal solution of equation (13.41) in the form of the series (13 .ip is unambiguously 
determined by the conditions of decoupling of correlations (I3.2l) - (l3.3p . Fourth, by the 
definition (I3.5p - (I3.6I) of functional C a {r, <p |R; z/} , the coefficients in its similar series 
expansion in terms of 0(r) satisfy the same conditions. To be exact, from (I3.2p - (I3.3I) 
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it follows that 



C n (r x ... r n | R; v) + y] / ... / C n+k (r 1 ... r n+k \ R; u) 

k=1 J n+1 J n+k 

if 9) (r| R; „) + £ ^ jf... jf C fc+1 (r, n ... r fe | R; „) 



at rj — > oo 



1 at R — >• oo , 



C(a, v) 

where 1 < j < n. At least, if we interpret the transitions (I3.2p - fl3.3l) in the sense of 
absolute integrability and speak about not too large density values (useful information 
in this respect can be found in [9]). 

Summarizing all the aforesaid, we can conclude that the solution to equation (13.71) 
is nothing but C CT {r, |R; is} = C{r, 0/(1 + er) | R; v(a, u)} . Combining this result 
with (I3.6P and (13.81) . we see that for arbitrary, in definite sense (see below), constant 
a and bounded function = 0(r) the equality 

z/C{r,a + 0|R; u} = C U | R; v(a,u) ) (3.9) 

takes place. It can be rewritten also as 

7>)C{r,0|R; v} = C (r, - 1 1 R; v(a,v)\ = C{r,0|R; v} 

(3.10) 

The left equality here defines the one-parametric family of such transformations of 
the functional C{r, | R; z/} which, according to the right-hand equality, conserve its 
value. One can easy verify that this family is the Abelian group with composition rules 

T(cr 2 ) T(cri) = T(cxi + cr 2 + a x a 2 ) , 
v(a 2 , v{ai , v)) = v(ai + a 2 + 0102 , v) (3.11) 

and with restrictions a > — 1 , 0(r) > —1 . The latter are clear in view of the fact 
that functional v [1 + 0(r)] C{r, | R; z/} represents mean concentration of atoms in 
presence of an external potential U (r) related to 0(r) via 0(r) = exp [ — U(r)/T ] — 1 
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(see e.g. [Hj). The substitution a = exp (a) — 1 yields T(a 2 )T(ai) = T(ai + a 2 ) 
and eliminates the restrictions. 

The infinitesimal equivalent of the equalities (13.91) or (13.101) looks most convenient 
when expressed via particular CF : 

|x(z/) + [l + x(z/)]z/^-|F 1 (e ' 7) (r|R;z/) = ujc 2 (r, r'|R; v) dr' , (3.12) 

[ 1 + h{v) ) v C n (ri ... r n | R; v) = u J C n+1 (r 1 ... r n , r'|R; z/) dr' , 



rax z/ 



X [V 



dC{a, v) 
da 



J <T=0 



^ / C 2 (r,0; z/)dr 



The function x(z/) defined just now is known (see e.g.[T2]) to be directly concerning 
the state equation of the system: 1 + h(v) = T(dv/dV)T , where V denotes the 
pressure. Notice that in the framework of the grand canonical ensemble very similar 
relations can be easy derived by differentiation of DF in respect to the activity. 



4. Invariance group of the generating functional 
of historical correlations 

We have approached to the substantiation of main results of the present work. Let 
us show that solution of the evolution equation (12.151) possesses invariance properties 
similar to (I3.10p . Since the initial condition to (12.151) (see (12.131) ) does not depend on 
the variables ip — ip(r, p) and v at all, the solution to (I2.15P is completely determined 
by the structure of operators C and C and the conditions at infinity (12.141) . The 
latter allow to extend the functional V{t, R, P, ip |Ro; v\ (like the functional C before) 
to arguments a + ip(r, p) , with a = const , in place of ip(r, p) . The fact that the limit 
in (I2.14p is achieved rapidly enough to indeed ensure this extension will become clear 
afterwards, from the relations (14.71) which will be obtained in this section. Next, thanks 
to (12.141) and besides to the definition (12. 6p of operator C , the variable ?/>(xi) in the 
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first term of the expression CV can be shifted by arbitrary constant: 

Pi d 5V f t N1 pi d 5V 



[ lj)(xx) — ■ 7^- ^ r = f[a + i){xx) ] — • 

Ji m or i dtp(xi) J x m ovi 



8ijj{x\ 



where a = const , for instance, a = 1 . That is important difference of CV from 
C T . Consequently, taking in mind the action of C onto V , one can write 

At last, let us look at the operator C ' (the formula next to (12.151) ). In contrast to £ , 
it depends on the density v . Nevertheless, with the help of equality (13.91) it is easy to 
make sure that it obeys a relation like (14.11) if the transformation of argument ip(r, p) 
is accompanied by transformation of argument v in accordance with (13. 8p : 

£ i"'° + *-£) = £ \^>Th'Wjh^]) (4 ' 2) 

The combination of all these observations implies the invariance property 

V{t, R, P, a + i> | R ; v] = V |t, R, P, | R ; v(a, v) | (4.3) 

or, equivalently and similarly to (13.101) . 

f(a) V{t, R,P, ^ | R ; u} = V{t, R, P, - 1 1 R ; v(a, v)} = 

1 + a 

= V{t,R,P, ^|R ; u} , (4.4) 

where the left equality, together with (13 .5h . ( 13. 8\) . ( I3.10P and (13. lip , defines action of 
the above characterized group of transformations onto functional V{t, R, P, ip \ R ; u} . 
Apparently, in respect to this functional that are such the transformations of its argu- 
ments which do not change its value. The expansion of the equality (14. 3p into series in 
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terms of ip yields 



K)(t,R,P|Ro; v(a,v)) = y (*,R,P|Ro; v) + 

+ E^- /•••/ K(t l R,rW | P l pW|H ;i/l 



(4.5) 



n = l 



n! 



1 



14(t,R,r( fc ),P,p( fc )|R ; ufai/)) = V k (t, R, r (fe) , P, p (fc) |R ; i/) + 



\4 +n (t,R,r( fc+ "),P,p( fc+ ™)|R ; v) (4.6) 



n = l ■-■ </fc+l J k+n 

Corresponding infinitesimal (in respect to a ) relations are similar to (13.121) 



(4.7) 



\nviy) + [l + x(z/)]z/|-| V n (t, R, r (n) , P, p (n) |R ; v) = 

= v [ K. + i(t,R,r (n+1) ,P,p (n+1) |R ; v) 

Jn+l 

Formulas ( 14.5ft and ( 14. 61) can be interpreted as "virial expansions" of the probabilistic 
law of BP's random wandering and statistical correlations between BP and medium, 
with those difference from the usual virial expansions of thermodynamic quantities p2] 
or kinetic coefficients [13] that here some decrements of the density do appear instead 
of its full value. However, in the limit when u — > and o — > oo with va = const 
our relations take quite usual form. For instance, ( 14.51) transforms into 

y (t,R,P|Ro;^)) = y o (t,R,P|Ro;0) + ^ ^ / K(t,R,r (n) ,P, P (n) |R ; 0) , (4.8) 



n = l 



V \V) = V 



n = l 



I/" 

n! 



1 Jn 



Cn+l( r ) r l ■■■ r n; 0) 



where J in the first row means integration over all the atoms' variables. It should 
be noted that formula (14. 5ft was obtained earlier in [TT1 [141 [T5] on the base of the 
"generalized fluctuation-dissipation relations" [TBI [17]. Up to that time, as far as I 
know, such kind of exact relations was not in use in statistical physics. It would be 
rather hard to extract them directly from the BBGKY hierarchy (12.21) or equivalent 
equations for CF. The only exclusion is the special case of BP in ideal gas [18], [T9]. when 
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C n — , and readers of this paper can easy derive relations (14 .7p (with x{y) = ) by 
means of differentiation of equations (12.171) in respect to gas density. 

5. On principal consequences from the invariance relations 

An exact solution of the BBGKY hierarchy automatically satisfies all the exact "virial 
expansions" (14 .51) - (14 M resulting from the above-stated invariance group. Moreover, 
expansion (14. 8 p itself gives explicit formal solution to the BBGKY hierarchy (since at 
v = the corresponding equations for CF can be solved by means of direct successive 
quadratures, like equations (I2.17P at v = ). Therefore it is natural to apply our above 
results to "testing of statistical hypotheses" about the BP's wandering. 

In this respect it is important to emphasize that formulas (I4.5l) -( j4~6l) interconnect 
two different random walks in two media with arbitrarily assigned densities ratio, 
v(a,v)/v. Consequently, all the terms in series (14 .5p or (14. 6p . as well as in (14. 8p . 
are equally necessary for correct summation of the series, and its truncation would 
produce not an approximative but wrong result. Even in the "dilute gas limit" (when 
gas parameters tend to zero, 47rrfz//3 — > , 47rr^z//3 — > with r a ^ being characteris- 
tic radii if atom-atom and BP-atom interactions) or in the so called "Boltzmann-Grad 
limit' (when simultaneously the mean free paths, = (irrlu)^ 1 and A a = {Kr 2 a v)~ x , 
are fixed), as well as in the case of ideal gas (when A a = oo but A;, is finite). In all 
these cases, formulas (14. 5ft - (14. 71) reduce to 

.J K, + ^R,r'-»,P,p<-»|R„;,),(5,) 

OV Jn+l Jn+k 

thus making it quite obvious that from the point of view of rigorous theory all the 
historical correlations are equally substantial. 

Meanwhile, it is the convention in physical literature to reject third- and higher-order 
inter-particle correlations when deriving closed equations (kinetic equations) for one- 
particle DF (for examples, see 0[TOl[T3],[2O]). Sometimes, as in the case of hard sphere 
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gas [20l[2T], these correlations get lost even insensibly for authors (see [18J about that). 
By the widely accepted opinion, under the Boltzmann-Grad limit the theory based on 
the BBGKY hierarchy should reduce to the classical kinetic theory of gases (where the 
BP's probability distribution, F Q = Vq , would obey either the linearized Boltzmann 
equation or similar Boltzmann-Lorentz equation [201 Ell [23] ) . However, really from the 
aforesaid it follows that such is not the case: the classical kinetics is not an extreme 
case of the statistical-mechanical theory. 

As an illustration, consider the BP in ideal gas. If rejecting all the higher-order 
correlations, i.e. putting on V2 = in the second of equations ( 12.171) (with n — 1 ) and 
then inserting its quadrature into the first equation, we come to a closed equation for 
Vq . An asymptotic of solution of this equation at t 3> r = Ab/vo (where vq ~ y/T/M 
is characteristic thermal velocity of BP) is certainly Gaussian: 



where AR = R — R , and D ~ t>oA{, is the BP's diffusivity. For a rarefied gas, 
A;, = {kt\v)~ x and D oc u^ 1 . Anyway, (15. 2p is some rather complicated function of 
v . At the same time, from the point of view of the exact relations (15.11) . the equality 
V2 = immediately implies that Vo(t, AR; v) is purely linear function of v ! 

So strong discrepancy does mean that asymptotic ( 15. 2ft is incorrect, and we have to 
get back to the complete BBGKY hierarchy. An approximate approach to its correct 
analysis was suggested already in [7] (or see [22]) and developed in [23]. In [23], for 
particular situation when the role of BP is played by a marked gas atom, the following 
asymptotic was found: 



V (t,AR;v)= / V (t,R,P\Ro;v)dP 



exp (—AR 2 /ADt) 
(AnDt) 3 / 2 



(5.2) 




(5.3) 
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Here, Q(x) is definite function which turns to unit at x = and rapidly goes to zero 
at x — >• oo . The characteristic feature of this probability distribution is that its fourth- 
order cumulants grow with time not linearly (as it would be in the case of Gaussian 
asymptotic) but under the law oc (Dt) 2 ln(t/r) , that is nearly proportionally to the 
square of time. 

Empirically, such the statistical property of the random walk is taken for the fluctu- 
ations in current diffusivity of BP [61 [21] whose power spectral density at frequencies 
/ <C 1/t is approximately inversely proportional to the frequency (so called "1/f -noise" 
[6]; thus, results of [3, [23] confirmed the conjectures about origin of 1/f -noise stated 
for the first time in [21] )• From the physical viewpoint, the origin of these fluctuations 
and corresponding "historical" statistical correlations described by CF V n ( n > ) is 
mere indifference of the system to the pre-history of BP's collisions, i.e. to their total 
number and their time-average rate, as well as to ratios of numbers of collisions with 
different impact parameters [7]. Thus, the historical correlations are caused by com- 
plicity of particles in their collision rate fluctuations rather than by their interactions 
as such. 

Let us demonstrate that quite similar conclusions can be deduced basing on only the 
virial expansions and the trivial non-negativeness of all the DF F n . 
From the non- negativity of F\ and identity (12.121) we have 

Vb(t,AR; v) [ Fi eq) (r\K;u)dr + [ K(t,R,r|Ro; v) dr > , (5.4) 
Jn Jn 

where Vi(t, R, r|R ; v) = J f Vi(t, R, r, P, p|R ; v) dprfP , and Q is any region in 
the space of vectors r — R . At given < 5 < 1 , let £1(5, t, AR; v) = Q(5) be the 
minimum (in the sense of volume) among all regions Q which satisfy the condition 



/ Vi dr - 




< 5 


J Vxdv 


Jn 
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From this definition it follows that 



Q max I V\ I > 



V x dr 



Vxdr 



(5.6) 



(spatial region and its volume are denoted by the same letter). In the aggregate the 
inequalities (15.41) and (15.51) produce, as can be easy verified, inequality 



V (t, AR; u) 



Q(S) 



Fi eq) (r\R;u)dr + (1 - 5) I \4(t,R,r|Ro; u) dr > 



(5.7) 



Combining it with the first of the virial expansions (14. 7\i (with n = ) and taking into 
account that in the case under consideration x(z/) — > and F^ eg \r\ R; v) < 1 , after 
quite obvious reasonings one obtains 

PiUJi A"R- 

(5.8) 



H(t,AR; ^)y (t,AR; ,/) + ^(t, AR; v) > () 



where the quantity f2(i, AR; z/) if defined as follows: 

n(5,t, AR; v) 



£l(t, AR; z/) = min 

0<5<1 



1-5 



(5.9) 



Next, discuss a meaning of the latter quantity. Notice, first, that Q(5) ^ when 
5 — > 1 , merely by the definition of Q(<5, t, AR; z/) = f2(<5) . Second, if the function 
Vi(t, R, r|R ; v) does not change its sign when changing argument r — R (but keeping 
other arguments fixed) then at 5 — > 1 the region Q(5) shrinks into infinitely small 
neighborhood of the point of maximum of |Vi| (or union of such neighborhoods), so 
that one can write 



1 V x dr - 1 V x dr 




J V x dr 


— Q(5) max Vx\ 


JQ(5) J 






r 



Q(S) max | Vi\ -> (1-5) 



Vxdr 



Consequently, Cl(t, AR; v) < \ J Vxdr |/(max r | VjJ) . From the other hand, because 
of the inequality definitely U(t, AR; v) > \ JVxdr\/( max r | Vx\) ■ From here 

(recalling the constancy of Vi's sign) we conclude that 



Q(t, AR; v) 





R 


\v)\dr 


max r 


^i(*,R,r 


Ro; u)\ 



(5.10) 
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This is nothing but the effective volume of a region in (r — R)-space on which the two- 
particle correlation spreads out. Now, let our assumption about the sign constancy of 
Vi(t, R, r|R ; v) in respect to r — R is not satisfied, and V denote those part of the 
(r — R)-space where V\ has dominating sign, that is coinciding with sign of integral 
J Vi dr . Then the analogous reasonings yield 
I / Vl dr 





JVidr 


^ Iv 


Vi\ 


dr 


max re x) 


Vi 


max 


r£T> 


\Vi\ 



. . < n(t,AR; i/) < |J ,' , < JPI ' (5.11) 
max r | Vi | " " 

Hence, in such the case we can say that Q(t, AR; v) is even less than effective volume 

of a region where significant two-particle correlation takes place. 

Further, let us confront inequality (15. 8p with the hypothesis that the asymptotic of 

Vo(t, AR; v) has Gaussian shape (15.21) . From (I5.8P it follows that this hypothesis can 

be true only if 

ATl , dh\D f AR 2 3\ 1 /AR 2 3\ /rioN 

n t.AR; u) > = (5.12) 

v ' ' ; - dv \ADt 2 J v \ADt 2 J v ; 

(we took into account that in gas D oc v~ l ). In other words, if the mentioned effective 

volume, occupied by the two-particle correlation, is not bounded above under variations 

of t and AR . In opposite, if this volume is bounded above, and 

i/0(t,AR; v) < C\ — const , (5.13) 

then the Gaussian asymptotic (15.21) is forbidden. Instead of it, inequality (15.81) allows 
the asymptotic like 

W' AR ^^(4^*(w) e (^) ' <5 ' 14) 



where function ^f(z) should satisfy the inequality 

z—^- + a^{z) > , a = - + ci 
dz 2 



dlnD 



-i -l 



3 
2 



ci (5.15) 



din v 

Consequently, ^(z) can not decrease at infinity faster than under the power law: 
ty(z — > oo) oc l/z a . The formula (15. 3p in [23] corresponds to c\ = uQ = 2 . 
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Thus, the theory inevitably leads to the statistical correlations whose length appears 
unbounded if not in space, as at variant (15.121) . then in time, as at variant (j5.13|) . 
It remains to ascertain what of these variants is more close to exact solution of the 
BBGKY hierarchy. From the physical viewpoint, the second is undoubtedly preferable. 
Indeed, according to equations (12.171) . the source of the correlations between BP and 
atoms is their collisions. A collision intends such relative disposition of BP and an 
atom when the vector p = r — R belongs to the "collision cylinder" which is oriented 
in parallel to the relative velocity u = p/m — P/M and has radius ~ . At the 
same time, the particles should not be separated by a distance much larger than A = 
min (A a ,Af,), since otherwise their collision most likely would be prevented by an 
encounter of one of them with other gas particles. Hence, the two-particle correlations 
are concentrated in the middle of the collision cylinder, at —A < u • p/\u\ < A . The 
volume of this space region does not depend on the particles' momenta, therefore it 
gives natural estimate of the volume Q(t, AR; u) . Obviously, it does not depend also 
on t and AR. Thus, Q ~ 2Ab7rr^ = 2jv (if, like on [23], BP is merely one of 
atoms). As the result, we come to (15.151) with c\ = uQ 2 . 

6. Conclusion 

To resume, the problem about thermodynamically equilibrium random walk of a 
test "Brownian" particle (BP) in a fluid was formulated as a problem of classical 
statistical mechanics in the framework of the Gibbs canonic statistical ensemble of 
initial conditions. The corresponding BBGKY equations for distribution functions 
(DF) and the Bogolyubov equation for their generating functional were considered 
in terms of the cumulant functions (CF) introduced in order to extract statistical 
correlations between current state of the fluid and total displacement, or path, of 
the BP during all the time of its observation ("historical correlations"). It is shown 
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that the generating functionals of equilibrium DF and time-dependent DF both possess 
invariance in respect to definite continuous group of transformations of their arguments, 
including the density of the fluid (mean concentration of atoms). 

The obtained invariance group produces the exact "virial expansions" ("virial rela- 
tions") which connect full sets of CF taken at different values of the density. With 
the help of these relations it was shown that all the many-particle correlations play 
equally important roles, regardless of the gas parameter value, and therefore the stan- 
dard "Botzmannian" approach to kinetics of BP, which rejects third- and higher-order 
correlations, is incorrect even in the low-density or Boltzmann-Grad limit. Correspond- 
ingly, the virial relations quite certainly forbid the Gaussian asymptotic of the BP's 
path probability density distribution inherent to the Boltzmann's kinetics. Instead, 
such an asymptotic is allowed only which has power-law long tails (cut off at the BP's 
path values on order of the ballistic flight path). This conclusion agrees with the re- 
sults of approximative solving of the BBGKY equations pO, [23]. It means that all the 
constituent parts of trajectories (random walks) of BP are statistically dependent one 
on another. Therefore, the trajectories cannot be imitated by sequences of mutually 
independent, in the sense of the probability theory, random trials (speaking figuratively 
and in the Einstein words, "God does not play dice" ). In this sense, any real molecular 
random walk, in contract to the random walk in the Lorentz gas, is non-ergodic. 

It should be noted that results of the present paper can be easily extended to ther- 
modynamically non-equilibrium random walks under influence of an external force (in 
this respect see [TU [15J, [22] ) . It would be interesting also to generalize the considered 
invariance group to another problems focused on not selected particles only but sooner 
on collective and hydrodynamical variables. 
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